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ABSTRACT: Two groups of parameters derived 
from both species diversity and species composition 
of spider assemblages of the East European Plain 
oak forests were tested in respect to large-scale 
variation in macroclimatic and habitat factors. The 
study is based on data equal in sample size and 
collecting efforts, Various diversity indices were 
taken as diversity measures. Spatial variation in the 
species composition was expressed by sample scores 
on the ordination axes derived from detrended 
correspondence analysis (DCA) and canonical 
correspondence analysis (CCA), both based on 
abundance and presence/absence data. Species 
composition is found to be a better indicator of 
spatial variation in macroclimate, giving higher and 
more stable values of correlation coefficients. 
‘Temperature parameters appear to be best explanatory 
variables for the patterns observed. The diversity 
measures vary with season and behave often in a 
different way for the spring, autumn or entire 
datasets. Only the total abundance shows a clear 
trend to a northerly and easterly decrease. 


РЕЗЮМЕ: Анализ зависимости двух групп 
параметров, вытекающих из видового 
разнообразия и видового состава группировок 
пауков в дубравах Русской равнины, от 
макрогеографической изменчивости климатичес- 


ких и биотопических факторов. Работа основана 
на данных, полученных по единой методике 
отбора проб и сбора материала. Для анализа 
биоразнообразия применялись различные 
индексы разнообразия. Пространственная измен- 
"HBOCTb видового состава оценивалась по 
характеру расположения учетов обилия H присут- 
ствия-остутствия на ординационных осях с помо- 
щью детрендного корреспондентного анализа 
(ОСА) и канонического корреспондентного 
анализа (ССА). Отмечено, что видовой состав 
лучше отражает пространственную изменчивость 
макроклиматических факторов, обнаруживая 
более высокие и стабильные показатели коэф- 
фициентов корреляции. Температурные парамет- 
ры оказались переменными, наилучшим образом 
объясняющие полученные данные. Меры био- 
разнообразия варьируют по сезонам, часто по 
разному отражая данные за весну, осень и весь 
полевой сезон. Лишь общее обилие обнаруживает 
четкую тенденцию к уменьшению к северу и 
востоку. 


Introduction 


The variety of diversity concepts has been 
divided by Whittaker [1977] into two main groups: 
(1) inventory diversity, and (2) differentiation 


! Contribution No, 3 to the project ‘Spatial variation in soil macrofauna communities of the East European oak forests in relation 
to environmental factors, sponsored by the former USSR Academy of Sciences 
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diversity. In his understanding, inventory diversity 
is intended to quantitatively describe the diversity 
of separate biotic assemblages (communities, 
ecosystems, biomes), while differentiation diversity 
assumes spatial aspects of variation in the diversity 
of biotic assemblages along gradients or among 
regions. In other words, inventory diversity describes 
the biotic assemblages irrespectively of the identities 
of the species they are composed of, meaning that 
two convergent communities situated in different 
continents may have equal values of various diversity 
statistics but completely different species listings. 
In contrast, differentiation diversity will depend on 
the similarity in the species composition of 
assemblages to be compared. The inventory diversity 
concept generally coincides with the landscape- 
typological approach as commonly termed in the 
Russian schools of ecology and biogeography, 
whereas differentiation diversity may be considered 
asan exampleofa regional-individualistic approach 
{see Chernov, 1975, 1984}. 

The above basic idea of Whittaker can be 
slightly modified and extended to cover any biotic 
parameters, be they strict diversity measures or not. 
Such an extension may be based on the borderline 
between the individualistic and typological 
approaches, that is, do the biotic characteristics 
consider the species composition or do they not. In 
this respect, two groups of parameters can be 
distinguished (сё. Hengeveld, 1990, Penev, 1992]: 

(1) Parameters dependent on the species 
composition expressed as species listings, based 
either on species presence/absence or abundance, 
frequency, density, biomass and/or other importance 
measures for each species. From these listings, a 
good number of structural characteristics can be 
derived, for example zoogeographical or phenological 
structures of communities. What is important, two 
assemblages situated in different faunistic or floristic 
regions will never be completely identical concerning 
such parameters. 

(2) Parameters independent of the species 
composition. They are derived irrespectively of 
species identities and include quantitative 
characteristics of the assemblage as a whole, such as 
inventory diversity statistics, total sum of individuals, 
total biomass, productivity, energy fluxes, etc. 
However, when based on functional classification 
principles or biological traits, structural 
characteristics are also to be grouped here, eg. 


trophic or life-form structures of communities. 

Usually papers dealing with spatial variation in 
diversity treat one of the two groups of parameters 
outlined above. This paper aims to compare species 
diversity with species composition in how both 
reflect large-scale variation of climate and habitat 
properties. This problem was discussed in an earlier 
paper based on soil-dwelling wire-worm assemblages 
(Coleoptera, Elateridae larvae) of oak forests along 
a latitudinal gradient in central Russia [Penev, 
1992]. Now we focus on spiders as a taxonomic 
group different from wire-worms in most biological 
properties, in particular as inhabiting an overground 
stratum, the soil litter. Two major gradients will be 
considered, one latitudinal reflecting trends in 
nature zonation from north to south, and the other 
longitudinal due to regional differences as well as to 
the increasing west-east climate continentality and 
severity. 


Material and methods 
Study area and field procedures 


Material was taken during five expeditions over 
the East European, or Russian, Plain carried out in 
May-June and from mid-August to early October 
1988 to 1990 (Fig. 1, Appendices 1-2). The oak 
forests were chosen as a model habitat since they 
represent the climax vegetation type in the central, 
nemoral regions of European Russia [Razumovsky, 
1981], thus being presumably the best indicators of 
macroclimate-determined spatial variation. The 46 
study sites, ranging over a distance of about 1,100 
km from north to south, are situated within three 
nature zones, from mixed coniferous-broadleaved 
forests to steppes. From west to east, the study area 
covered a distance of about 2,300 km, ranging from 
Moldova (= Moldavia) to the Urals. 

So many study sites scattered over such vast 
distances can hardly be expected to be too similar 
to each other. Indeed, as shown in Appendices 1-2, 
the oak forests chosen occupy different relief forms, 
they are supported by a variety of soil types, and 
they can be represented by various plant associations. 
Besides, there is no doubt that local habitat 
conditions do influence large-scale variation in 
biotic assemblages, modifying the major role of 
climate. To minimize between-site variation, we 
sticked to the most mature and well-preserved 


Fig. 1. Geographical location of the study sites. Site numbers as in Appendices 1-3. Solid lines show borders between zones, as 
follows: TAI — taiga, МХ — mixed coniferous-broadleaved forest zone, FST — forest-steppes, STE — steppe zone. 

Puc. 1. Географическое местоположение точек сбора материала. Номера точек как в Приложениях 1-3. Сплошные линии 
показывают границы между зонами: TAI - тайги, MIX - смешанных хвойно-широколиственных лесов, FST - лесостепи, STE - степи. 
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forest patches in each site, representing plant 
associations most typical for that particular area. In 
central Russia, for example, zonal oak forests on 
interfluve plains are largely Quercetum aegopodiosum 
and Quercetum aegopodioso-caricosum pilosae with 
their variants, in southern Moldavia these are 
Quercetum pubescentis cotinosum and Quercetum 
pubescentis herbosum; in Cisuralian lowlands they 
are mainly Quercetum mixtoherboso-aegopodiosum 
[Nikolaeva, 1963; Geideman, 1968, Gorchakovsky, 
1972, Kleopov, 1990]. Apart from this, the sites 
were chosen so as to avoid marginal variants as 
regards moisture conditions, that is both too wet or 
too dry oak forests were not included in the field 
work. 

Twenty soil samples of 25 x 25 cm each were 
taken at each site. The soil animals were collected 
from 10 cm thick soil layers on oilskin in the field 
according to Ghilarov [1975]. The spiders were 
present only in the litter and were sampled by 
means of a sieve. The samples were taken along 
short transects (200-500) m within more or less 
homogeneous forest patches in consequent series of 
3 samples each: the first sample at an oak trunk, the 
second between trees, and the third in a small forest 
‘window’. Other details with regard to sampling 
sites/methods can be derived from Penev [1992]. 


The environmental data set 


There are two major climatic gradients affecting 
the distribution of oak forests on the Russian Plain. 
The first one can be called north-southern and it 
reflects latitudinal changes in thermal conditions 
and in the radiation balance as a whole [Grigoryev, 
1954, Budyko, 1977]. The second gradient is 
orientated from west to east, being largely due to 
the increased climatic continentality and severity. 
To trace these two major trends, we chose climatic 
parameters (Appendix 3) correlated either with 
latitude (TANN, TMIN, FRL, PR, WA, CO HUMI, 
DEF) or with longitude (TAMP). Most temperature 
parameters, however, decrease not only from south 
to north but also in a west-east direction. 

Some of the environmental variables were 
transformed according to the following scales: 

1) Mean annual precipitation (PR) (mm): 1 — 
400-500, 2 — 500-600, 3 — 600-700. 

2) Precipitation for the warm period (WA): 1 — 
250-300, 2 — 300-350, 3 — 350-400, 4 — 400-450, 
5 — 450-500, 6 — 500-550. 

3) Precipitation for the cold period (CO): 1 — 
100-125, 2 — 125-150, 3 — 150-175, 4 — 175-200, 
5 — 200-225. 


4) Mechanical composition of the soil (MECH): 
1 — light sandy, sandy loam or loess-like soils, 2 — 
moderate loam, 3 — heavy loam to clayish; in 
borderline cases 0.5 is added to the lower value. 

5) Thickness of the litter layer (LIT) (cm): 1 — 
0-2, 2 — 2-4, 3 — 4-6, 4 — 6-8. The values in 
Appendix 3 were estimated as an average of the 20 
soil samples taken at each site. 

6) Humus depth (HUM) (cmy: 1 — 0-5, 2 — 5- 
10, 3 — 10-15, 4 — 15-20, 5 — 20-25, 6 — 25-30, 7 
— 30-35, 8 — 35-40, 9 — >40. The values in 
Appendix 3 represent the averages for each site 
estimated from the 20 soil samples. 


The species data set 


The species composition and abundances at 
each site are given in Appendices 4-5. Out of a total 
of 452 spider species known from oak forests of the 
East European Plain in general [Esjunin et al., 
1993], 153 species were found in our soil samples. 
All calculations, however, were carried out 
considering only adults representing but 124 species, 
all typical inhabitants of the soil litter stratum. 
Since we found considerable variation in the 
samples taken in the spring and in the autumn 
[Esjunin et al, 1993], the statistical procedures 
were applied to the whole data set, as well as to 
subsets including either autumn or spring sites. 


Data analysis 
Species composition 


The sample scores on the first four DCA axes 
computed with CANOCO [Ter Braak, 1988] were 
considered as a measure of between-site variation in 
the species composition. For the same purposes, we 
also applied a direct gradient analysis technique, i.e. 
Canonical Correspondence Analysis, ССА [Ter Braak, 
1985]. Because of a strong multicollinearity among 
most of the climatic factors, a forward selection of 
the environmental variables was carried out before 
the computation of CCA axes. We chose the 
variables with the highest values of variance explained 
andstopped when the Monte Carlo permutation test 
(99 random permutations) showed no significant 
additional effect of the variable in question. To 
better reflect the existing environmental dataset, we 
also selected some soil variables with rather low 
values of explained variance. The significance of the 
correlations between DCA axes and each 
environmental variable were estimated in all cases by 
means of transformation of the correlation coefficient 
to Fisher's z-statistics [Lakin, 1990]. 
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Fig. 2. Single linkage (nearest neighbour) clustering dendrogram of Kendall's rank correlation coefficient between quantitative 
characteristics of the spider assemblages: S — number of species, N — individuals total, Ц — Melntosh's diversity measure and D — 


Melntosh's D index (Magutran, 1988, р. 40), Dmg — Margalef's index, Dmn — Menchi 
index (index of ‘polydominance’) (= Hill's N,), Sh = exponential Shannon's function ( 


index, 1/C — reciprocal of Simpson's 
= Hill's М), H' — Shannon's index, PIE — 


‘probability of interspecific encounter’, NMS — number of moves per specimen, HB — Brillouin's diversity index, HBE — Brillouin's 
зев, E — Pielou's evenness measure, E" — Hill's evenness measure (= М), Е — Alatalo’s evenness, Ө — Molinari's evenness, 
C — Simpson's index (‘concentration of dominance’), d — the Berger-Parker index (= reciprocal of Hill's №). 

Puc. 2, Кластерная дендротрамма одноразовой связи (ближайшето соседства) коэффициента рантовой корреллуни Kenaasaa 
для количественных характеристик группировок пауков: $ - число видов, М - суммарное число экземпляров, D - коэффициент D 
Макинтоша [Magurran, 1938: 40], Dmg - индекс Маргалефа, Dan - индекс Менхиника, 1/C - величина, обратная коэффициенту 
Симпсона (индекса “полиломинантности") (= индекс N, Хиллз), Sh - экспонента коэффициента Шеннона (= индекс N, Xanna), 


W' - коэффициент Шеннона, PIE - “вероятность межвидовой nerpem”, NMS = число перемещения на экземпляр, НВ - индекс 


разнообразия Брийлузка, HBE - индекс выровененности Брийлузна, 


Е - индекс выровененности Пайлоу, E” - индекс выровсненности 


Хилла (= Мун, F = индекс выровененности Алатало, G - индекс выромеменности Молинари, C - коэффициент Симпсона 
(“концентрация доминирования”, 4 - индекс Бергера-Паркера (= величина, обратная индексу N, Хилла). 


Species diversity 


The choice of an appropriate index of diversity 
is a common and not yet fully solved problem in 
modern ecology. An incredible amount of work, 
based on both theoretical and empirical evidence, 
have been devoted to this problem fi.e. Hutcheson, 
1970; Hurlbert, 1971; Hill, 1973; De Benedictis, 
1973; Peet, 1974; Pielou, 1975; Routledge, 1979; 


Alatalo, 1980; Pesenko, 1982; Pesenko & Semkin, 
1989; Magurran, 1988; Molinari, 1989, etc.]. In any 
event, we rather prefer to stick to 'pluralistic' views 
(eg. Hurlbert, 1971; Pesenko, 1982; Pesenko & 
Semkin, 1989] and use several indices, although 
most of them are known to more or less strongly 
correlate with each other. Our choice of diversity 
indices was based on the following presumptions: 
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(1) Theoretical properties of the indices 
depending on what aspect of diversity they are 
designed to measure; 

(2) Empirical conclusions drawn from a 
correlation analysis of the various indices based on 
our field data. 

We shall avoid giving here all the formulas of 
diversity indices, since they are easily available 
from several recent reviews fe.g. Pesenko, 1982, 
Magurran, 1988]. Besides, their abbreviations are 
already commonly used in the modern ecological 
literature. 

As shown by Hill [1973], many diversity indices 
are members of the equation: 


N-(p*p,*p,*-p)s 


where а > 0. When а = 0, №, = 5 (number of species); 
when a = 1, М, = Sh (exponent of the Shannon 
function; when a = 2, №, = 1/С (inverse of 
Simpson's 'concentration of dominance’ index), etc. 
Thus, with an increase in а, the dependence on the 
number of species is decreased. The last member of 
Hill's series, the reciprocal of the Berger-Parker 
index 4, is completely independent of 5 [ Magurran, 
19887. Hence, we chose both the first three and the 
last member of Hill's equation. In addition, some 
other diversity statistics were tested, namely 
Brillouin’s function (ИВ), Menchinic’s (Dmn) and 
Margalefs (Dmg) indices. 

Evenness measures have extensively been 
discussed in several papers [Hill, 1973; Alatalo & 
Alatalo, 1980; Alatalo, 1981; Molinari, 1989]. Hill's 
evenness E’ (= N,/N,) has been criticized by Peet 
[1974] and Molinari [1989], since it may have high 
values in cases when the real evenness is low. 
Instead, the measure F has been proposed [Alatalo 
& Alatalo, 1980; Alatalo, 1981]. However, in its 
turn, Alatalo’s index was criticized by Molinari 
[1989] who proposed a calibrated variant of this 
index, called G measure. Based on the above 
considerations, we selected the indices F, G, as well 
as Brillouin's evenness HBE. 

The dominance measures usually negatively 
correlate with diversity. However, their use 
sometimes allows to illustrate cases of high 
predominance of one or two species. The most 
popular dominance measure is Simpson's index C. 
The Berger-Parker index can also be considered as 
a dominance measure, since it merely expresses the 
relative abundance of the most abundant species. 

The correlation structure of the diversity indices 
applied is shown in Fig. 2. The nearest-neighbour 
clustering dendrogram based on Kendall’s rank 
correlation coefficient was used as an additional 


criterion for omitting some highly correlated 
diversity statistics such as the ‘number of moves per 
specimen’, NMS [Fager, 1972; Pesenko, 1982], the 
‘probability of interspecific encounter, PIE, 
Shannon's and McIntosh’s diversity indices, H’ and 
U respectively, Pielou’s and Hill's evenness indices, 
Е and Е” respectively. 

The calculation of all diversity statistics and the 
cluster analysis were performed with the program 
BIODIV [Baev & Penev, 1991]. 


Results 


Environmental factors underlying variation in 
Species composition 


The DCA ordination results are presented in 
Figs 3a-f, the correlations of the DCA axes with 
each environmental variable in Tab. 1, the 
eigenvalues and percentage variances accounted for 
by the first two axes are given in Tab. 2. 

The ordination scattergrams based on the 27 
spring siles (Figs 3a, 3b) demonstrate the major 
role of macroclimatic gradients in compositional 
variation of the spider assemblages. The distribution 
of siles on Axis 1 for presence/absence data 
illustrates an isolated position of the Cisuralian oak 
forests (sites 42, 44, 46), on the one hand, and of 
those lying in the central part of the mixed 
coniferous-broadleaved forest zone (sites 18, 19, 
22), on the other. In general, the first DCA axis can 
hardly be considered asa clear geographical gradient, 
although it correlates with geographical longitude. 
Axis 1 shows a relationship to temperature, while 
Axis 2 to humidity parameters (Tab. 1). The 
scattergram based on abundance data has a similar 
structure and correlates with the same factors (Fig. 
3b, Tab. 1). 

The dominant compositional gradient for the 
autumn datasct on Axis 1 is orientated latitudinally 
and correlates positively with the mean annual 
temperature and the deficit in air humidity (Fig, Зс, 
Tab, 1). The second axis, in its turn, shows a 
significant correlation with precipitation parameters. 
For the abundance data (Fig. 3d), Axis 1 arranges 
the sites from southwest to northeast and strongly 
correlates with latitude, longitude, and all 
temperature parameters as well as with the 
mechanical composition of the soil. The second axis 
showed a complicated gradient which correlating 
with none of the environmental variables. 

The first DCA axis generated on the basis of the 
entire (spring + autumn) dataset can be interpreted 
as a clear latitudinal gradient for both presence/ 
absence (Fig. 3e) and abundance (Fig. 3f) data. The 
values of the correlation coefficients (Tab. 1) 
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Table 1. Correlations between the assemblage parameters and environmental factors, Abbreviations of environmental 
variables as in Appendix 2; ~ — not significant, the numbers with an asterisk (°) are significant at P < 0.01, those 
with two ones (“*) at P < 0.001, the remaining at P < 0.05. The assemblage parameters are abbreviated as follows: C 
— Simpson’s index (‘concentration of dominance’), D — McIntosh’s D index, d — the Berger-Parker index (=Hill’s 
N), DCA1(a) — first DCA axis based on abundance data, ОСА1 (р) — first DCA axis based on presence/absence 
data, DCA2 — second DCA axes respectively, Dmg — MargalePs index, Dmn — Menchinic’s index, Е — Alatalo's 
evenness, G — Molinari’s evenness, HB — Brillouin’s diversity index, HBE — Brillouin's evenness, N — individuals 
total, 5 — number of species, Sh — exponent of the Shannon's function (= Hill's N,), 1/С — inverse of Simpson's 
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Таблица 1. Корреляции между параметрами группировок и факторами среды. Сокращения переменных 
окружающей среды как в Приложении 2: ‘-’— недостоверно; числа, снабженные звездочкой (^), достоверны при 
P < 0.01, числа, снабженные двумя звездочками (7), достоверны при P < 0.001, прочие при P < 0.05 Параметры 
группировок сокращены: С — коэффициент Симпсона (“концентрация доминирования”), D — коэффициент D 
Макинтоша, d — индекс Бергера-Паркера (= индекс М, Хилла), ОСА1(а) — первая ось детрендного 
корреспондентного анализа (AKA) на основе данных по численности, ПСА1 (р) первая ось ДКА на основе 
данных по присутствию/отсутствию, ОСА? — соответственно вторые оси ДКА, Dmg — индекс Маргалефа, Dmn 
— индекс Менхиника, F — индекс выровененности Алатало, G — индекс выровененности Молинари, НВ — индекс 
разнообразия Брийлуэна, HBE — индекс выровененности Брийлуэна, М — суммарное число экземпляров, $ — число 
видов, ЗВ — экспонента коэффициента Шеннона (= индекс №, Хилла), 1/С — величина, обратная коэффициенту 
Симпсона (индекса “полиломинантности”) (= индекс М, Хилла) 
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Table 2. 


Eigenvalues and percentage variance accounted for by the first two DCA and CCA axes. 


Таблица 2 


Собственные значения и процент изменчивости за счет первых двух осей детрендного 


и канонического корреспондентного анализа 


Eigenvalue Percentage variance Eigenvalue Percentage variance 
DCA1 DCA2 DCA1 DCA2 ССА1 CCA2 CCAT CCA2 
Spring samples 
Presence/absence 0.57 035 85 3.8 0.49 0.37 209 153 
Abundance 0.75 0.36 122 69 0.64 0.37 28.0 159 
Autumn samples 
Presence/absence 0.53 042 12.3 97 0.48 0.40 209 173 
Abundance 0.51 0.39 125 8.5 0.47 0.37 21.5 170 
All samples 
Presence /absence 0.54 0.33 13.4 61 0.41 0.32 20.2 158 
Abundance 0.50 0.39 110 73 0.44 0.37 228 18.9 


suggest that climate, ie. primarily temperature, 
precipitation and air humidity, can be considered as 
the major factors affecting the large-scale variation 
in spider assemblages. It is noteworthy that 
significant correlations of the site scores with the 
thickness of the litter stratum for both presence/ 
absence and abundance data are also observed 
(Tab. 1). The arrangement of sites on both second 
axes reflects a longitudinal variation in climate, 
expressed by an increasing temperature amplitude 
and decreasing remaining temperature parameters 
(presence/absence data). The mechanical 
composition of the soil and its correspondent 
variable ‘sandy and sandy loam soil’ correlated with 
Axis 2 for abundance data (Tab. 1), suggesting that 
soil properties may influence species abundances, 
perhaps indirectly via the plant cover and/or the 
composition of forest litter. 

In most cases, except for the autumn dataset 
(see Tab. 1), we did not find significant correlations 
in the third and fourth DCA axes which would add 
much to the understanding of the observed patterns. 
The third axis for abundance and, partly, presence/ 
absence data correlated with soil factors, such as the 
‘mechanical composition of the soil’, ‘presence of 
sand’, ‘podzolized soil type’, and ‘humus depth’. 

In general, the analysis based on all three 
datasets showed that temperature parameters are 
the most important determinants of the main trend 
in the species data, regardless of its latitudinal or 
longitudinal orientation. Among the temperature 
measures, the mean annual temperature and the 


number of frostless days seem to have the most 
stable and highest values of the correlation 
coefficient. The humidity parameters can be judged 
as being of subordinate importance. The soil 
properties seem to influence the species abundances 
by modifying the large-scale variation caused by 
climate. In both entire and autumn datasets, the 
species composition based on presence/absence 
appears to be largely influenced by macro- 
environmental factors, whereas the species 
composition combined with abundances reflects 
local soil conditions as well. 

Direct gradient analysis (Figs 4a-f, Table 3) 
confirms the major role of temperature. In four of 
the six datasets, the variable with most important 
contributions to the multiple regression (i.e. with 
the highest T-value of the canonical coefficient on 
the first CCA axis) is either the number of frostless 
days or the mean annual temperature. The most 
important determinants for the spring dataset are 
the number of frostless days, the mean annual 
temperature amplitude, and the thickness of the 
litter stratum (Figs 4a, b). The relative air humidity 
correlates with the second CCA axis. The autumn 
samples (Figs 4c-d) show quite a similar pattern 
demonstrating the importance of the mean annual 
temperature (or its exchangeable variable, the 
number of frostless days), the annual temperature 
amplitude, the relative air humidity (or its 
counterpart, the deficit in air humidity). The role 
of podzolized soil type on both biplots (Figs 4c-d) 
can be considered as related to the geographical 
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Fig. 4. CCA biplots of sites and environmental variables. Designations as in Fig. 3, variables as in Appendix 2. 


Рис. 4. Биплоты KaHORMNCCKOTO корреспондентного анализа точек сбора материала и переменных окружающей среды. 
Обозначения как на Рис. 3, переменные как в Приложении 2. 
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Table 3. Canonical and correlation coefficient of the CCA axcs and environmental variables remained after 
CANOCO forward selection. The values of the canonical coefficient marked with an asterisk (^) have T-values 
greater than 2.1 for both spring and whole datascts and 2.23 for autumns samples respectively. The variables not 
selected for the corresponding dataset are designated with "-" 

Таблица 3. Канонический коэффициент и коэффициент корреляции осей канонического корреспондентного 
анализа и переменных окружающей среды, оставшиеся после отбора программой CANOCO, Всличины 
канонического коэффициента, отмеченные звездочкой (“), имеют значения Т свыше 2,1 как для весенних 
учетов, так и дая всего банка данных, и 2,23 дая осенних учетов. Переменные, не отобранные для 
соответствующего банка данных, обозначены как “-” 


Г мене — — — | — ner ins 
Canonical coefficients Correlation Smee ile 
en qr x 
Spring samples 
FRL -0.56* 0.27* -0.86 -001 -1.07* 022 -0.81 -015 
|TAMP 0.22* 0.37* 0.68 0.29  -011 0.53* 0.49 018 
HUMI 0.15* 0.46* 024 0.59  -007  101* 0.04 0.87 
co -0.08 -0.04 -0.16 0.02 -0.24* -0.01 -0.22 0.01 
LIT 0.25* 0.17 -0.50 -0.04 014 035* -0.49  -037 
HUM -0.39* 0.23* -0.27 043  -047* -005 -0.27 0.17 
MECH 0.17* -0.05 0.11 0.02 0.18 0.15 0.13 0.23 
PODZ 0.07 0.26+ 0.08 041 0.03 0.30% -0.03 0.32 
CHER 0.01 0.01 -0.11 0.11 -0.01 -0.08  -014 -0.24 
Autumn samples 
FRL -0.12 -0.29* 0.59 0.13 - = - - 
ITANN - - - - -0.65* -0.03 -0.96 0.05 
TAMP -0.47* -0.25* -0.16 -0.68 0.08 -0.05 0.70 0.46 
DEFI 0.72* 071* 0.81 -0.26 - - - - 
HUMI - - - - -0.10 -045* 004  -0.87 
WA -0.18 0.79* -0.40 0.72 0.20 0.08 -0.01 -0.65 
LIT -0.12* 0.14 0.12 0.17 0.05 -0.01 -0.07 0.05 
HUM -0.25* 0.10 -0.01 -0.27 005 -0.26* 0.01 -0.17 
MECH -0.01 0.31* -0.28 -0.03 -0.03 -0.05 0.47 0.15 
PODZ -0.07 0.19 -0.30 0.56 -0.07 -036* 014 -0.68 
All samples 
TMIN -0.02 0.28 0.60 0.64 - - - - 
TANN - - - - -0.74* -0.03  -090 002 
TAMP -0.26 -0.04 -0.30 -0.77  -048 0.06 0.63 0.56 
DEFI 0.65* 0.02 0.88 -0.22 - - >. - 
HUMI - - - - 0.08  -048* 017  -090 
WA -0.06 0.36* -0.51 0.66 = - - - 
co - - - -0.13* 0.04 0.02 -049 
LIT -0.09 0.12* 0.52 0.08 -0.02 -0.01 -0.38 030 
HUM -0.06 0.11* 0.16 -0.02 -0.17* -0.05 -0.12 017 
MECH -0.01 0.17* -0.20 -0.04 002 0.04 0.36 0.25 
PODZ 0.02 0.16* -0.36 0.50 -019* -0.12* -0.05 -066 
GREY 0.03 -0.04 0.09 -0.38 - - > - 
CHER - - - - 0.05 0.13* -0.07 042 
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trend in the species data, since this type of soil is 
characteristic of the northern sites. The whole 
dataset demonstrates much of the same pattern 
(Figs 4e, Гу. Precipitation for the warm period 
having high importance in Figs 4c and e, we consider 
this as a result of the influence of a single deviating 
site (Fig. 1: site 3), the Rostochye Reserve, situated 
close to the foothills of the Carpathians and 
characterised by high precipitation rates. 

The length of the arrows on the CCA biplots 
suggests that the species composition expressed in 
a presence/absence form shows a stronger relation 
to climatic parameters in comparison to species 
abundances. The latter seem to be influenced not 
only by climatic variables but also by factors 
characterizing the properties of the soil. It is 
rather possible that some noise in the species 
abundance data can also be attributed to seasonal 
variation. 


Patterns of variation in species diversity 


Similar to the species composition, the species 
diversity parameters were tested for correlation 
with environmental variables separately for the 
three datasets (Tab. 1). 

The spring samples demonstrate an increase in 
the total abundance and diversity (number of 
species 5, indices 1/C, Sh, and HB) both from north 
to south and from west to cast. This pattern is 
illustrated by significant negative correlations of 
the above parameters with latitude, longitude, 
annual temperature amplitude as well as with a 
decrease in the remaining temperature factors. On 
the contrary, the dominance indices show a reverse 
pattern meaning that, along with an increase in the 
climate continentality and severity, the 
‘concentration of dominance’ within the spider 
assemblages grows, while their diversity becomes 
lower. Evenness measured with the indices Fand G 
seems to be poorly related to macro-environmental 
factors and positively correlates solely with the 
number of frostless days. 

The total abundance of the autumn samples 
correlates in the same way with climate, i.e. it is 
decreased both castwards and northwards. The 
diversity measures, however, show a reverse pattern 
in comparison to the spring samples (Figs 5a-b), 
that is they are increased with latitude and longitude. 
Similarly, the evenness grows northeasterly, being 
negatively correlated with temperature and positively 
with the annual temperature amplitude. On the 
contrary, the more southern spider assemblages 
demonstrate a higher ‘concentration of dominance’ 
in the autumn. 


It is noteworthy that the spider assemblages on 
soils of a heavier mechanical composition have a 
lower total abundance and a higher evenness as 
judged from autumn samples. However, this pattern 
is not observed in the spring dataset. In both 
autumn and spring samples, both the total abundance 
and the diversity as measured with Brillouin's index 
correlate positively with the thickness of the litter, 
although Brillouin's index in itself is known to be 
sometimes related to the total abundance [Pect, 
1974]. Significant correlations of the evenness with 
the humus depth is revealed for both spring and 
whole datasets. 

Analysis of the entire dataset confirmed strong 
negative correlation of the total abundance with an 
increase in the climate continentality and severity. 
Diversity measured with Brillouin's function also 
correlates positively with temperature. In particular, 
some of the evenness measures (HBE, F) as well as 
the diversity index D demonstrate a weak but 
significant relationship with both the mechanical 
composition of the soil and the humus depth, 
suggesting that the richer and the heavier the soil, 
the more diverse the assemblages. Similarly, the 
sites with a thicker litter show higher diversity as 
measured with the index ИВ. 

The discriminant abilities of two popular indices, 
the № measure of diversity and the E measure of 
evenness (Pielou’s evenness index), are illustrated 
in Figs 6a-b. It is evident that the evenness measure 
fails to reveal any patterns of geographical variation 
in the assemblages as compared to № (see also Fig. 
5). The same conclusion can also be drawn for the 
№/№ measure of evenness and the Shannon function, 
none indicating a dependence of the spider 
assemblages on the environmental variables. 


Discussion 


Species diversity versus species composition 
as a reflection of large-scale variation in climate 
and habitats 


The above results allow to conclude that the 
relation of diversity statistics to environmental 
variation is seasonal. Moreover, although being 
more or less correlated to each other, most of the 
diversity measures react differently according to 
the environmental variables. In the spring and 
autumn, the diversity is decreased and increased 
eastwards, respectively. Evenness fails to demonstrate 
a clear geographical trend in the spring, but in the 
autumn it is likewise increased to the cast. The 
differences among the indices cannot be considered 
as pure mathematical artifacts, they may reflect real 
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Fig. 5. Scattergrams of the М, diversity measure (= reciprocal of Simpson's index) of the spider assemblages. The location of the 
sites corresponds to the DCA scattergrams in Figs ЗЬ and 2d, respectively. 


Puc. 5. Скаттерграммы индекса разнообразия N, (= величина, обратная коэффициенту Симпсона) группировок пауков. 
Расположение точек соответствует скаттерграммам ДКА соответственно на Puc. 3b и 2d. 
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Fig, 6. Discriminant abilities of the N, diversity measure (= reciprocal of Simpson's index) and the E evenness measure (= Pielou's 
index of evenness) based on our filed data. The location of the sites corresponds to the DCA scattergrams in Fig. ЗЕ 


Рис. 6. Дискриминаитные возможности индекса разнообразия N, (= величина, обратная коэффициенту Симпсона) и индекса 


выровненности E (коэффициента выровненностя Пайлоу) на основе наших даиных. Расположение точек соответствует 
скаттерграммам AKA на Рис. 3f. 
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Fig. 5. Seattergrams of the N, diversity measure (= reciprocal of Simpson's index) of the spider assemblages. The location of the 
sites corresponds to the DCA scattergrams in Figs 3b and 2d, respectively, 


Рие, 5. Скаттерграммы индекса разнообразия N, (= величина, обратная коэффициенту Симпсона) группировок пауков. 
Расположение точек соответствует скаттерграммам AKA соответственно на Рис. ЗЬ и 24. 
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Fig. 6. Discriminant abilities of the N, diversity measure (= reciprocal of Simpson’s index) and the E evenness measure (= Pielou’s 
index of evenness) based on cur filed data. The location of the sites corresponds to the DCA scattergrams in Fig. 3f. 


Рис. 6. Дискриминантные возможности индекса разнообразия М, (= величина, обратная коәффициенту Симисона) и индекса 


выровненности E (коэффициента выровнеиности Пайлоу) на основе наших данных. Расположение точек соответствует 
скаттерграммам ДКА на Рис. ЗЕ 
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underlying processes in spider assemblages. For 
example, the spring samples from the southwestern 
oak forests are characterised by the highest species 
richness. In the Moldavian sites (10, 11, 13, 14), 14- 
21 species per site were found; in the sites of the 
Dnieper basin (4, 6, 7, 8) 11-26 species per site; in 
the central Russian oak forests (18, 19, 20) 7-16 
species per site; and in the Cisuralian ones (42, 44, 
46) 8-12 species per site. The high species richness 
of the southwestern oak forests is accompanied by 
codominance of 2-3 or even 4 species: Microneta 
viaria, Diplocephalus picinus, rarerly Enoplognatha 
ovata (juveniles), Tapinocyba insecta, Diplostyla 
concolor, Lepthyphantes flavipes, Hahnia nava, Nigma 
flavescens. At the same time, solely juveniles of 
Helophora insignis (sites 44, 46) or Ceratinella brevis 
(site 42) dominated the Cisuralian oak forests. 

In the autumn, the southwestern oak forests are 
characterised by a high-level species richness. The 
diversity, however, is decreased at the expense of 
the increasing predominance of a single species, 
Microneta viaria, while in the sites located within 
the zone of mixed coniferous-broadleaved forests 
(16, 17, 21, 23, 24) no prevailing single-species 
dominance is observed. In the Cisuralian sites (40, 
45), the dominance of Helophora insignis revealed in 
the spring is smoothed in the autumn by its 
codominance with other species such as Lepthyphantes 
menget, L. tenebricola, Bolyphantes alticeps. 

The above considerations suggest that species 
diversity can be regarded as a discrete community 
parameter influenced by a larger number of stochastic 
factors than species composition. Variation in 
diversity among local assemblages depends on 
climate, habitat peculiarities, species phenology, 
possible biotic interactions and historical reasons. 
Although the diversity statistics, primarily the total 
abundance and the number of species, display some 
trends in geographical variation, the patterns 
observed are more irregular in comparison to 
species composition. Moreover, it is evident that 
the conclusions made on the basis of one or a few 
diversity statistics should be interpreted but with 
caution. Hence, the diversity indices seem to be less 
suitable as indicators of macro-environmental 
changes, particularly of climate. This conclusion is 
further reinforced by evidence recently presented 
by Helióvaara et al. [1991]: "Species richness 
decreased northwards, but the relation was not very 
convincing, suggesting that more specific factors 
may affect diversity in the Scolytidae" (p. 17). 

On the contrary, compositional changes along 
vast geographical gradients seem to have much in 
common with the differentiation of biotic 
assemblages on a local scale. А unimodal response 


of a species to environmental factors [Whittaker, 
Ter Braak, 1990], usually proved for local habitat 
gradients, may be extended on a far larger 
geographical scale. Hengeveld [1985] found that 
the Dutch carabid beetles are normally distributed 
in relation to climatic variables. Later he proposed 
the species-surface model [Hengeveld, 1990] which 
predicts that the abundance of a certain species will 
gradually decrease from the centre of the range to 
its margins. Therefore, local species assemblages 
will be composed of species deriving from different 
parts of their ranges. Chernov [1989] supposed that 
continualistic views on spatial variation of a biota 
should be extended to a zonal scale. Hence, species 
responses to climate on broader scales will result in 
continual variation patterns in the assemblage 
composition, similar to the individualistic concept. 
of Ramensky-Gleason created to explain plant 
species distributions along local gradients. On the 
other hand, the climatic parameters also vary 
continually and therefore should be better indicated 
by continually varying biotic parameters, such as 
species composition, rather than discrete ones, such 
as species diversity. However, the above conclusion 
presumes that this approach may be applied only 
within more or less homogeneous faunistic/floristic 
regions. For example, a continual latitudinal variation 
in species assemblages along the entire north-south 
gradient from the European tundra up to the alpine 
belt of the Caucasus can hardly be expected with 
respect to climate, sincea highly significant faunistic 
boundary has long been revealed between the 
Caucasus and Europe at which the fauna drastically 
changes due to merely historical reasons. In such 
cases, the second group of parameters (independent 
of the species composition) provides the only basis 
for comparisons in assemblages responses to the 
recent environmental conditions. 


Environmental determinants in large-scale 
variation in species assemblages 


The importance of temperature for the 
distribution of animals has long been known (e.g. 
Hesse, 1924; Uvarov, 1931]. Over the last few 
decades, an increased interest has appeared to 
broad-scale effects of thermal conditions on the 
distribution of biotic assemblages as well. Latitudinal 
variation in species richness has been discussed in 
numerous papers [reviews: Brown & Gibson, 1983; 
Hengeveld, 1990], being usually attributed to 
global changes in the radiation balance from the 
equator to the poles. 

It has been shown by many authors that the 
diversity parameters of spider assemblages on local 
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scales correlate with the thickness and structure of 
forest litter [Lowrie, 1948; Huhta, 1971; Uetz, 
1976, 1979; Bultman et al., 1982; Leclerc & Blandin, 
1990]. Our results demonstrate this relationship on 
a far broader geographical scale, too. 


Conclusions 


This study aimed to prove one important 
methodological problem, namely which species 
assemblage parameters better reílect large-scale 
variation in climate. Besides, some patterns in 
spider assemblage variation with respect to the 
environment were also discussed. The following 
conclusions can be drawn as regards both issues: 

(1) In studies on the geographic variation in 
spider assemblages, seasonal changes within them 
should be accounted for. This results from biological 
peculiarities of spiders, particularly from their 
adults often restricted to relatively short time 
periods. The same conclusion can perhaps be 
inferred for all soil invertebrate groups which 
support no more or less constant abundance levels 
of their populations during the whole year, either 
separately for imaginal and larval stages or both. 

(2) The major variation trend in the species 
composition can be attributed to temperature, 
regardless of a latitudinal or a longitudinal 
arrangement of sampling sites. This conclusion 
holds for the autumn, spring, and whole datasets as 
well as for both presence/absence and abundance 
data. Among the temperature measures, the mean 
annual temperature and the number of frostless 
days show the highest values of correlation 
coefficients. The humidity parameters (precipitation, 
relative air humidity and its deficit) as well as soil 
properties (mechanical composition, thickness of 
litter, humus depth) seem to be of subordinate 
importance, 

(3) North-south trends in the faunal composition 
can be related to climate, while those from west to 
east may be considered as a result of the combined 
effects of both recent ecological conditions 
(increasing climate continentality and severity) 
and historical processes (glaciations, influence of 
two important nemoral refuges, the Carpathian and 
the Uralian - see Esjunin et al, 1993). 

(4) Species richness is the highest in the 
southwestern oak forests, decreasing both northerly 
and easterly. This pattern, however, cannot be 
regarded as a regular gradient and is not supported 
by stable significant correlations with environmental 
variables. 

(5) The total abundance of the assemblages 
decreases significantly north- and eastwards with 


the increasing climate continentality and severity. 

(6) Numerous alpha-diversity indices also vary 
in relation to climatic and habitat factors, but in 
different, sometimes contrary, ways in the spring 
and in the autumn. 

(7) The parameters derived from species 
composition seem to be better indicators of spatial 
variation in climate rather than the parameters 
independent from species identities, e.g species 
diversity. Species composition varies continually 
within relatively homogeneous faunistic regions, 
thereby reflecting continual variation in climate. 
Species diversity can be regarded as a discrete and 
stochastic biotic parameter strongly influenced by 
a large number of factors. 

(8) Species composition based on presence/ 
absence data largely reflects macroclimatic gradients, 
while species listings combined with abundances 
demonstrate also a relationship with local habitat 
conditions, primarily with the mechanical soil 
composition and litter thickness. 
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№ 
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te 


Struki 12 km NE Buda-Koshelevo 


Derno 30 km E Lutsk 


Rostochye Reserve ca. 40 km 
WNW Lvov 


Peruny пг. Mikhaylovka 
10 km $ Kanev 


Peruny nr. Mikhaylovka 
10 km $ Kanev 


Chernyi Les nr. Zuamenka 
Komissarovski Les 12 km W Volnogorsk 
Savranski Les nr. Polyanetsko 
Savranski Les nr. Polyanetsko 
Rossoshanskii Les 3 kin N Briceani 
Kodry Reserve nr. Lozova 
Kotovskii Les nr. Kotovsk 
Kotovskii Les nr. Kotovsk 
Tigechskii Les 5 km $ Tigecei 
Sokolovo 15 km E Dubniki 
Dubniki 16 km NINE Udomlya 
Sharapovo 50 km SW Moscow 
Sharapovo 50 km SW Moscow 
Vasilevskii 10 km N Veney 
Tulskiye Zaseki nr. Krapivna 
Tulskiye Zaseki nr. Krapivna 
Chernysheno ca. 30 km WSW Kozelsk 
Trud са. 20 km SE Ulyanovo 
Setukha 50 km E Orel 

Khinel 25 km SW Sevsk 

Kazatskii Les 25 km SW Kursk 


Les-na-Vorskle Reserve 


Province Date 
Gomel 7-9.09.1989 
Volyn 12-14.09.1989 
Lvov 16-20.09.1989 

Cherkassy 28-29.05.1988 

Cherkasy 28-29.09.1990 

Kirovograd 20-22.05.1988 


Dnepropetrovsk 19-21.05.1990 


Odessa — 31.05-2.06.1988 
Odessa 4-7.10.1989 
Moldova 29-31.05.1990 
Moldova 7-9.06.1988 
Moldova 26-29.09.1989 
Moldova = 1-206.1990 
Moldova 5-6.06.1990 
Novgorod 31.08-3.09.1989 

ee 4-6.09.1990 
Moscow 9-10091990 
Moscow 23-24.05.1991 
Tala 7-10.05 1988 
Tula 3-7.05.1989 
Tul 12-13,09.1990 
Каша 2804-10589 
Kaluga — 18-2208.1991 
Orel 15-16.09 1990 
Bryansk — 18-19.09.1990 
Kursk 5-9.05.1990 
Belgorod 15-20.05.1988 


Appendix 1 
Geographical location and characteristics of 46 study sites (see also Fig. 1). When not stated otherwise, the 
dominant oak species is Quercus robur |. Soil properties in Appendix 2. 
Приложение 1. 
Географическое местоположение и характеристика 46 точек сбора (см. также Рис. 1). Tac не указано иначе, 

ломинирующей породой дуба выступает вид Quercus robur L. Свойства почв приведены в Приложении 2. 


Forest type 


interfluve Querceto-Carpinetum 
nudum 


interflave Querceto-Carpinetum 
nudum 


upland Quercetum nudum 


Quercetum aegopodiosum on fluvial 
terrace 
Quercetum aegopodiosum on fluvial 
terrace 
interfluve Fraxincto-Quercetum 
aegopodiosum 
interfluve Querceto-Fraxinetum 
mixtoherboso-urticosum in ravine 
interfluve Quercetum 
mixtoherboso-aegopodiosum 
interfluve Quercetum 
mixtoherboso-aegopodiosum 
interfluve Quercctum 
mixtoherboso-urticosum 
upland Quercetum petraeae 
aegopodiosum 
upland Quercetum pubescentis 
cotinosum 
upland Quercetum pubescentis 
cotinosum 
upland Quercetum nudum 


interfluve Quercetum 
mixtoherboso-oxalidosum 


foodland Quercetum caricosum 


interfluve Quercetum 
aegopodioso-caricosum pilosae 
interfluve Quercetum 
aegopodioso-caricosum pilosae 
interfluve Querceto-Tilietum 
equisetosum 
interfluve Querceto-Tilietum 
aegopodiosum 
interfluve Querceto-Tilietum 
aegopodiosum 
interfluve Quercetum 
aegopodioso-alliosam 
interflave Quercetum 
aegopodioso-galeobdolosum 
interfluve Quercetum mixtoherbosum 
in ravine 
Quercetum caricoso-aegopodiosum on 
fluvial terrace 
interfluve Quercetum 
aegopodioso-convallariosam 
upland Quercetum 
aegopodioso-caricosum pilosae 


120-130 


100-150 
80-120 
80-120 

100-120 
100-200 
80-120 

120-200 

120-200 
70-120 
80-100 

40-70 
40-70 

70-150 

200-220 

150-250 

120-150 

120-150 
50-70 

200-300 

200-300 

-100150 

150-400 

100-120 

100-170 
50-100 


200-300 
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Appendix 1 
Geographical location and characteristics of 46 study sites (see also Fig. 1). When not stated otherwise, the 
dominant oak species is Quercus robur L. Soil properties in Appendix 2. 


Приложение 1. 
Географическое местоположение и характеристика 46 точек сбора (см, также Рис. 1). ac He указано иначе, 


доминирующей породой дуба выступает вид Quercus robur L. Свойства почв приведены в Приложении 2. 


Forest 
No te Province Date Forest type 
age (уг) 
1 Struki 12 km NE Bada-Kosheleva Gomel 7-9.09.1989 Бене ове тшен 120.130 
2 Derno 30 km E Lutsk Volyn 12-14.09.1989  interfluve аасы 100-150 
3 ate esee ca. 40 km Lvov 16-20.09.1989 upland Quercetum nudum 80-120 
Peruny nr. Mikhaylovka. | Quercetum aegopodiosum on fluvial 
A saint Chee Cherkassy — 28-2905.1988 AS 80-120 
5 Perany nr. Mikhaylovka Cherkassy — 2829091999 Quercetum acgopodiosum on avial уду 
10 lan $ Kanev zm d terrace 
6 Chernyi Les nr. Znamenka Kirovograd 20-22.05.1988 "erfluve Fraxineto-Quereetum әр 
acgopodiosum 
7 Komissarovski Les 12 km W Voluogorsk Dnepropetrovsk 19-21.05.1990  terfuve Querceto-Fraxinetum o 10 
mixtoherboso-urticosum in ravine 
$ Savranski Les nr. Polyanetsko Odessa 31.05-2.06.1988 interfluve Quercetum 120-200 
mixtoherboso-aegopodiosum 
9 Savranski Les nr. Polyanetske Odessa 47.10.1989 interfluve Quercetum 120-200 
mixtoherboso-aegopodiosum 
= М interfluve Quercetum " 
10 Rossoshanskii Les 3 km N Briceani Moldova 29-31.05.1990 Grebe eens 70-120 
11 Kodry Reserve nr. Lozova Moldova — 79.06.1988 upland Quercetum petraeae 80-100 
acgopodiosam 
12 Kotovskii Les nr. Kotovsk Moldova — 26-29.09.1989 "Pland Quercetum pubescentis 40-70 
cotinesum 
13 Kotovskii Les nr. Kotovsk Moldova 1-2.06.1990 upland Quercetum pubescentis 40-70 
cotinosum 
14 Tigechskii Les 5 km $ Tigecci Moldova 5.606990 apland Quercetum nudum 70-150 
а interfluve Quercetum 
15 Sokolovo 15 km E Dubniki Novgorod — 31.08-3.09 1989 анана 200-220 
16 Dubniki 16 km NNE Udoml „уе 4-6.09.1990 floodland Quercetum сайсозаш 150-250 
ya (=Kalinin) 
interfluve Quercetum E 
17 Sharapovo 50 km SW Moscow Moscow 9-10091990 aam pose 120150 
E m А р interfluve Quercetum 
18 Sharapovo 50 km SW Moscow Moscow 23-24051991 О em pilose 120-150 
19 Vasilevskii 10 km М Venev Tula 7-1005.1988 а 50-70 
equisetosum 
reeto-Til 
20 Tulskiye Zaseki nr. Krapivna Tula 3-7.05.1989 EM 200-300 
“Тен 
21 Tulskiye Zaseki nr. Кгарува Tula 12-13.09.1990 ме ат к 200-300 
es 2 2 interfluve Quercetum 
22 Chernysheno ca. 30 km WSW Kozelsk Kaluga 280410389 ры -100150 
р РЕ " Е 5 interfluve Quercetum 
23 Trud са. 20 km SE Ulyanovo Kaluga — 18-2208.1991 зен женге А 150-400 
24 Setukha 50 km E Orel Orel 15216091990 interfluve Quercetum mistoherboraam 100-120 
25 Khinel 25 km SW Sevsk Bryansk 18-19.09.1990  Queteetum у» трон ?^ 100.170 
A ER interfluve Quercetum 5, 
26 Kazatskii Les 25 km SW Kursk Kursk 5-9.05.1990 АЕ зет е ыйгар 50-100 
27 Les-na-Vorskle Reserve Belgorod 15-20.05.1988 upland Quercetum, 200-300 
aegopodioso-caricosum pilosae 


———  scgopocioso-carkosum мо | 
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Appendix 1. (Continued) 
Приложение 1. (Продолжение) 


№ Site Province Date Forest type pet 
age (уг) 
28 Les-na-Vorskle Reserve Belgorod — 24-2609.1990 upland Quercetum 100-300 
acgopodioso-caricosum pilosae 
29 10 km SW Trostyanets Sumy 21-23.09.1990 — iuterfluve Quercetum aegopodiesum 100-150 
30 Venevitinovo 33 km NE Voronezh Voronezh — 2930051989  Querecto-Pinetum acgopodioso- 70-80 
caricosum (pilosae) on fluvial terrace 
31 Voronezh Reserve Voronezh 31.05-1,06,1989 _ Quereetum acgopodioso-caricosam сд 
pilosae on fluvial terrace 
32 Voronezh Reserve Voronezh 23.06.1989 Quercetum aegopodioso-caricosum 70_100 
pilosae on fluvial terrace 
33 Tellerman Forestry nr. Borisoglebsk Voronezh 20-24.05.1989 interfluve Quercetam 100-150 
aegopodioso-caricosum pilosae 
34 Shipov Les nr. Vorontsovka Voronezh  5-6.06.1989 upland Quercetum 80-130 
acgopodioso-caricosum рйозас 
35 Shipov Les nr. Vorontsovka Voronezh 7-8.06.1989 upland Quercetum 80-150 


aegopodioso-caricosum pilosae 
36 Gaidary $ Gotvald Kharkov 12-13.05.1990 — interfluve Quercetum aegopodiosum 120-150 
floodland Quercetum 


A PEE K _ - 

37 Serebryanskoye Forestry nr Kremonnaya Lugansk — 1517051990 MALE 80-180 

38 Stanichno-Laganski Reserve nr. Lugansk Lugansk 11-13.06.1989 foodland Quercetam 70-90 
convallarioso-aristolochiosum 

39 Volga-Kama Reserve са. 50 km S Kazan Tartar 9-1 floodiand Quercetum 70-140 


mixtoherboso-acgopodiosum 


40 Zhigali Reserve, nz. Bakhilovo Samara 23-26.08.1988 __ Upland Quereetum 80-100 
mixtoherboso-aegopodiosum 


41 Movosemeikino ca. 15 km NW Samara Samara 18-2208.1988  interfuve Quercetum convallariosum 60-80 

42 Sarashi 20 km $ of Barda Perm 25-26.06.1991 olan Керен 60-120 
aegopodioso-graminosum 

43 Sarashi 20 km $ of Barda Perm 4-509.1991 upland Quercetum 60-120 
acgopodioso-graminosum 

44 Vilyai 20 km $ Asha Chelyabinsk 15.06.1991 _ interflave Quercetum 60-80 

mixtoherboso-aegopodiosum 
"mon vm upland Quercetum Р 
45 Shulgan-Tash Reserve Bashkiria 12-15.08.1988 Кена едда 100-200 
46 Shulgan-Tash Reserve Bashkiria 17.06.1991 upland Quereetum 100-200 


mixtoherboso-aegopodiosum 
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Appendix 2. Soil properties (sce Appendix 1). 
Приложение 2. Свойства почвы (см. Приложение 1). 


Soddy podzolic 
2 Soddy podzolic Sandy loam 15-30 
3 Soddy podzolic to carbonatic Sandy and loess-like loam to. gravel 8-25 
4 Soddy podzolic Sandy loam to sandy 3-12 
5 Soddy podzolic Sandy loam to sandy 12-30 
6 Grey forest Moderate loam 10-20 
7 Chernozem Loess-like to moderate loam 30-40* 
8 Grey forest Moderate loam 15-35 
9 Grey forest Moderate loam 8-35 
10 Grey forest Moderate loam 30-40* 
11 Brown forest Moderate loam 10-30 
12 Grey forest Loess-like loam 5-15 
13 Grey forest Loess-like loam 3-12 
14 Brown forest Sandy & loess to heavy loam 2-20 
15 Soddy podzol gleyed Clayish 5-15 
16 Soddy podzolic gleyed Heavy loam 10-20 
17 Soddy podzol Moderate loam 8-12 
18 Soddy podzol Moderate loam 10-25 
19 Grey forest Moderate to heavy loam 20-30 
20 Grey forest Moderate to heavy loam 10-25 
21 Grey forest Moderate to heavy loam 10-25 
22 Grey forest, partly gleyed Moderate to heavy loam 10-20 
23 Grey forest Heavy loam 5-15 
24 Chernozem Moderate loam 10-50* 
25 Grey forest Sandy loam 15-20 
26 Chernozer: Moderate loam 20-40* 
27 Grey forest Loess-like to heavy loam 10-40* 
28 Grey forest Loess-like to moderate loam 10-40* 
29 Grey forest Moderate loam 15-30* 
EU Grey forest podzolized Sandy loam 3-20 
31 Grey forest Moderate to sandy loam 8-40* 
32 Grey forest Silty to sandy loam 5-40* 
33 Grey forest Moderate to heavy loam 25-40* 
34 Grey forest Moderate to heavy loam 22-40* 
35 Grey forest Moderate to heavy loam 7-40* 
36 Grey forest Moderate loam 8-30* 
37 Meadow chernozem-like Moderate to heavy loam 20-40* 
38 Meadow chernozem-like podzolized Loess-like to moderate loam 10-40* 
39 Soddy podzolic Moderate to sandy loam 8-30 
40 Brown forest Moderate loam, gravel 5-40* 
41 Chernozem Moderate loam 15-40* 
42 Grey forest Gravel 3-20* 
43 Grey forest Gravel 1-20 
44 Grey forest Moderate loam, gravel 5-20 
45 Grey forest Moderate loam, gravel 10-40 
46 Grey forest Moderate loam, gravel 7-20 
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Appendix 5. (Continued) 
Приложение 5. (Окончание) 


1 2 » 5 9 12 15 1é 17 2 23 24 25 28 29 40 41 43. 45 
$ 14 15 21 11 17 17 24 15 11 9 7 14 18 21 20 re 15 17 16 


N 19 163 99 60 101 117 060 41 20 21 24 21 40 130 80 58 35 27 55 
D 05 057 069 05 07 065 082 079 081 08 064 09 083 058 074 075 O81 085 077 
Dmg 272 275 435 244 347 336 562 3.77 3.34 263 189 427 461 411 434 394 394 485 3.74 
Dmn 128 117 211 142 1.69 157 31 234 246 196 143 306 285 184 224 223 254 327 216 
1/C 414 451 684 31 734 596 125 909 741 723 411 1131 1081 445 849 829 935 9.99 914 


Sh 626 646 1027 $06 1033 871 1724 1134 8.98 7.93 5.23 12.59 1381 8.41 1204 11.88 11.54 13.32 11.48 


НВЕ 0.65 064 07 052 074 07 083 08 079 09 076 088 087 068 083 082 086 084 082 
0.6 0.64 063 052 068 064 071 078 08 09 074 089 077 047 068 067 079 073 078 
021 027 025 014 031 027 035 045 048 064 039 062 045 01 O31 03 046 038 044 


032 0.14 017 008 O11 O14 014 024 009 0.09 0.22 012 012 011 Q1 011 
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